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ABSTRACT 

We apply population synthesis techniques to calculate the present day number of two 
types of white dwarf-main sequence star (WDMS) binaries within the cataclysmic 
variable 2 to 3 hour period gap. The first are post-common envelope binaries with 
secondary stars that have masses 0.17 < AI2/MQ < 0.36 (gPCEBs), such that they 
will commence mass transfer within the period gap. The second type are systems that 
were CVs at some point in their past, but detached once they evolved down in orbital 
period to « 3 h as a consequence of disrupted magnetic braking, and are crossing the 
period gap via gravitational radiation (dCVs). Full population synthesis calculations 
are performed where we assume either constant, global values of the common envelope 
ejection efficiency, acE, or consider acE as a function of secondary mass. Several forms 
of magnetic braking are also considered. We predict an excess of dCVs over gPCEBs 
within the period gap of ~4 to ~13 assuming acE = 0.1 — 0.6, and an initial mass ratio 
distribution of the form n{q) = 1. This excess is revealed as a prominent peak at the 
location of the period gap in the orbital period distribution of the combined gPCEB 
and dCV population. We suggest that if such a feature is observed in the orbital period 
distribution of an observed sample of short orbital period WDMS binaries, this would 
strongly corroborate the disruption of magnetic braking. 
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1 INTRODUCTION 

The standard paradigm of cataclysmic variable (CV) evolu- 
tion assumes that orbital angular momentum (AM) is lost 
via both gravitational radiation and magnetic braking at 
orbital periods Porb/h > 3, and purely by gravitational ra- 
diation for Po,b/h<3 (e.g. Warner 1995; King 1988). How- 
ever, the AM loss rate associated with magnetic braking is 
ill-constrained, and it is unclear if magnetic braking ceases 
to be a sink of AM for Porb/h <3 (tKolb,2002i). 

Extrapolating the work of lSkumanichI (|l972l ') who inves- 
tigated the spin dow n rates of G-type stars by a magnetically 
coupled stellar wind, IVerbunt fc ZwaanI (|l98ll ) obtained the 
magnetic braking AM loss rate of the form Jmb oc Q,"^ , where 
n is the spin frequency of the donor star. 

As a consequence of the AM losses and the mass loss 
they drive, the donor star is unable to maintain thermal 
equilibrium. This causes donors wi th masses < 0.8 Mp, 
to ex ceed their equilibrium radius l|Stehle. Ritter fc Kolbl 
Il996l ). Mass loss will eventually cause the donor star to be- 
come fully convective at a mass M2 = M2,conv and when 
Porb /h ~ 3. At this point it is suggested that either the mag- 



netic activity of the donor star discontinuou sly decreases, or 
the topology of the magnetic field changes (|Taam fc SpruitI 
1989). This is based on the result that the interface between 
the radiative core and the convective e nvelope play a key role 
in th e generation of the magnetic field (|McDermott fc TaamI 
Il989h . 

When magnetic braking ceases at M2 = M2,conv the 
donor shrinks within its Roche lobe, mass transfer is shut 
off and the secondary re-attains thermal equilibrium. The 
system is no longer observed as a CV. Subsequent orbital 
evolution of the CV is driven by gravitational radiation only. 
Mass transfer resumes again when the secondary star makes 
contact with its Roche lobe once more at Porb/h ~ 2. 

This disrupted magnetic model was conceived to ex- 
plain the lack of observed CVs in the range 2 < Porb /h ^ 3 
(Rappaport, Verbunt & Joss 1983; Spruit fc Ritter 1983), as 
shown in Figure [1] 

However, despite the success of this model in repro- 
ducing the correct width and position of the period gap 
(e.g. Kolb, King fc Ritter 1998), there is no independent 
observational support for disrupted magnetic braking. Ob- 
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Figure 1. Orbital period distribution of CVs from Ritter & Kolb 
(2003), catalogue edition 7.7 (2006). The axis on the left indicates 
the number of systems per bin. The solid, red curve is the cumula- 
tive distribution with numbers, TV < Porbi indicated on the right 



servations of young, open clusters (Sills, Pinsonneault & 
Terndrup 2000) show that the spin down rate of rapidly 
rotating single stars is much less than predicted by the 
IVerbunt fc Zw aan^ (!l98ll) relation. Also, no discontinuous 
decline in magnetic activity for the fully convective stars 
is detected (Jones, Fischer & Stauffer 1996). It is likely that 
for the high rotational speeds of > 100 km s~^ of the CV 
donor stars, the AM loss rate, Jmb, due to magnetic braking 
will saturate (f or a review see e.g. Collier-Cameron 2002). 
ISills et al] l|200Qi ) suggest the form Jmb (x in the saturated 
regime. Within this scheme the AM loss rate in CVs would 
be a factor of 100 less than that for the standard magnetic 
braking model, and the donor star would not be sufficiently 
drive n out of thermal equilibrium and no p eriod gap would 
arise l|Andronov. Pinsonneault fc SiUs|[200a ). 

Such theoretical and observational problems with the 
standard model of CV evolution aside, it is the only one 
thus far that predicts a period gap of the correct width with 
edges at the correct orbital periods Q. Two recent pieces 
of evidence in support for the disrupted magnetic braking 
model lies in the mass-radius relation of donor stars in CVs 
(Smith & Dhillon 1998; Patterson et al. 2005; Knigge 2006) 
and in the effective temperatures of the accretion-heated 
white dwarfs in CVs (Townsley & Bildsten 2003). In the 
former, the donor mass-radius relation for Porh /h ^ 3 lies 
systematically above that for isolated main sequence stars. 
This is consistent with the donor star being driven out of 



1 Willems et al. (2005) and Willems et al. (2007) studied CV 
populations evolving under the influence of circumbinary disks 
and found a period gap feature, but had to appeal to additional 
mechanisms in the pre-CV evolution to account for the lower edge 
of the gap. 



thermal equilibrium. Furthermore, a discontinuous decrease 
in the donor radii between 2 < Porb /h ^ 3 has been detected. 
In the latter piece of evidence, the effective temperatures 
of white dwarfs are higher above the period gap than be- 
low, and are overall consistent with the expected accretion 
heating for mass transfer rates predicted by the disrupted 
magnetic braking model. 

Politano fc Weiler (2006) calculated the secondary star 
mass distribution of present day PCEBs for a range of 
AM loss prescriptions, including disrupted magnetic brak- 
ing. They found that the number of PCEBs per secondary 
mass interval decreases by 38 per cent once magnetic brak- 
ing begins to drive their evolution for secondaries with ra- 
diative cores. They suggest that this decrease in the num- 
ber of PCEBs with secondaries between ~ 0.25 and 0.5 M© 
provides an observable test for the disruption of magnetic 
braking. 

In this paper we suggest an alternative, even more di- 
rect, test for the disruption of magnetic braking by investi- 
gating a more explicit consequence: the existence of a pop- 
ulation of WDMS binaries that were CVs in the past but 
became detached at Porb/h « 3 and are currently crossing 
the period gap. We call these dCVs. 

Using population synthesis techniques we calculate the 
present day population of these dCVs, as well as WDMS 
binaries that have emerged from a common envelope (CE) 
phase. We call these systems post-CE binaries (PCEBs). 
In particular, we only consider PCEBs which will begin 
mass transfer within the period gap. For our stellar models 
and choice of gap boundaries, such systems have secondary 
masses in the range 0.17 < M2/M-Q < 0.36. This sub-set 
of PCEBs we henceforth designate as gPCEB systems. Note 
that gPCEBs have Porb/h > 2 but can also have Porb/h > 3. 

As we will show below orbital the period distribution 
of the combined population of dCVs and gPCEBs reveals 
a prominent peak within the period gap due to an ex- 
cess of dCVs over gPCEBs here. We calculate the ratio 
dCV:gPCEB within the period gap, which indicates the 
height of this peak, for a range of models concerning the val- 
ues of the CE ejection efficiency, acE, the initial distribution 
of secondary mass, and the AM loss rate due to magnetic 
braking. The detection of such a feature in the orbital pe- 
riod distribution of an observed sample of WDMS binaries 
with orbital periods of a few hours would provide a simple 
test for the disruption of magnetic braking. 

The Sloan Digital Sky Survey (SDSS) is currently cre- 
ating the opportunity for large-scale observational studies 
of WDMS binaries, having identified already more than 
1000 systems (Silvestri et al. 2007; Schreiber, Nebot Gomez- 
Moran & Schwope 2007). Attempts to establish the subsam- 
ple of PCEBs among those systems and measuring their or- 
bital periods as well as their stellar properties is underway 
(Rebassa-Mansergas et al. 2007, Schreiber et al. 2008). We 
are optimistic that a quantitative test of the calculations 
carried out in this paper, and hence of the validity of the 
disrupted magnetic braking hypothesis, will be feasible in 
the near future. 

The structure of the paper is as follows. In Section 2 we 
describe the computational method, in Section 3 we present 
our simulation results, which are discussed in Section 4. We 
conclude our investigation in Section 5. 
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2 COMPUTATIONAL METHOD 

The computational method we apply consists of two main 
steps; the first is to calculate the zero-age populations of the 
dCVs (i.e. CVs that have just detached at the upper edge of 
the period gap) and PCEBs (i.e. WDMS binaries that have 
just emerged from a CE phase) using the population syn- 
thesis code BiSEPS (WiUems & Kolb 2002; Willems & Kolb 
2004). These populations are unweighted i.e. no formation 
probabilities are yet calculated for the systems. The code 
e mploys the single star evolut ion (SSE) formulae described 
in iHurlev. Pols fc ToutI l|2000l ) and a binary system evolu- 
tion s cheme based on that described in lHurlev. Pols fc ToutI 
l|2002h . The second step uses a different code introduced 
in Willems et al. (2005), which calculates the present day 
populations of these systems using a library of evolutionary 
tracks, here produced by BiSEPS itself. The procedure is 
described in more detail in the following sections. 



2.1 Initial Binary Population 

BiSEPS evolves a large number of binary systems, which 
initially consist of two zero-age main sequence stellar com- 
ponents. The stars are assumed to have a population I chem- 
ical composition and the orbits are circular at all times. The 
initial primary and secondary masses are in the range 0.1 
to 20 Mq, while the initial orbital periods range from 0.1 to 
100 000 d. There is one representative binary configuration 
per grid cell within a three-dimensional grid consisting of 
60 logarithmically spaced points in primary and secondary 
mass and 300 logarithmically spaced points in orbital pe- 
riod. There are ~ 5.4 x 10® binaries that are evolved for a 
maximum evolution time of 10 Gyr. For symmetry reasons 
only systems with Mi > M2 are evolved. 

We calculate the four dimensional probability density 
functions -^wd, Af2, Porb) of all zero-age dCVs and 
(j}{t, MwD, M2, Poih) for all zero-age PCEBs using the ini- 
tial distribution functions quoted below. Here t is the time 
since the birth of the Galaxy, Mwd is the mass of the white 
dwarf, M2 is the mass of the secondary star and Porb is 
the orbital period. The number of zero-age main sequence 
primaries, dA*', that form with masses in the range dA/i is 
given by dA'^ — /(Mi)dMi, where the probability density 
function, /(A/i), is given by 





^^^^'^ - ^ 0.15571M,-'-2 
0.15571M-^-'' 



A'/i/Mq < 0.1, 
0.1 < Mi/Mq < 0.5, 
0.5 < Mi/Mq < 1.0, 
1.0 < Mi/Mq, 



(1) 



(Kroupa, Tout & Gilmore 1993). The probability that a zero- 
age main sequence secondary star forms with a mass A'h 
can be determined from the initial mass ratio distribution 
(IMRD) given by 



n(g) = 







< g < 1, 
q > 1, 



(2) 



where q = M2/M1, is a constant, and fi is a normalisation 
factor such that J n{q) dq — 1 {/i = u + 1 for u > —1). We 
also consider the case where the secondary mass is deter- 
mined from an IMF according to equation ([T} where Mi is 
replaced with M^. We use n{q) = 1, i.e. u = 0, as our ref- 
erence model. Finally, the probability that a binary forms 



with an initial orbital separation a is determined by (Iben 
& Tutukov 1984; Hurley et al. 2002) 



h{a) = 





0.078636a- 



a/RQ < 3 or u/Rq > 10^ , - 
3<a/RQ<10^ ^' 



We further assume that all stars in the Galaxy are 
formed in binaries and that the Galaxy has an age of 10 
Gyr. The star formation rate, 5", is calculated by assum- 
ing that one binary per year is formed with Afi > 0.8 M0. 
Combining this with the Galactic volume of 5 x 10^^ pc'^ 
gives an average local birthrate of white dwarfs of 2 x 10~^^ 
pc~^ yr"'^, consistent with observations (Weidemann 1990). 
The lower limit of 0.8 Mq is the smallest mass the primary 
star can have if it is to evolve into a white dwarf within the 
lifetime of the Galaxy. We therefore have 



f{Mi)dMi = lyr" 



(4) 



From ([1} and Q we obtain S = 7.6 yr~^. We assume that 
this is constant throughout the lifetime of the Galaxy. 

We then employ a population calculator to convert 
x(t, A/wD, A^2,Porb) and (f>{t, Mwd, M2, Porb) into distribu- 
tions over Mwd, M2 and Porb describing the present day 
population of dCVs and gPCEBs. The secular evolution of 
each system in the zero-age populations is therefore followed 
using BiSEPS, and its contribution to the population, deter- 
mined by the above distribution function and star formation 
rate, is accounted for at each stage of its evolution. 



2.2 Binary Evolution 

The CE phase (|Paczvnskilll976r ) is modelled by equating the 
binding energy of the primary's envelope to the change in 
the total orbital energy of the binary system. Quantitatively 
this is 



GM,{M, + A/c 

APl.L 



acE 



GM0M2 G{Mc + M,)M2 



2at 



2ai 



(5) 



where G is the gravitational constant. Ale and Mc axe the 
masses of the primary's core (the proto-white dwarf) and 
envelope, Pi,l is the radius of the primary immediately be- 
fore the onset of the CE phase, and M2 is the mass of the 
secondary. The orbital separation of the stellar components 
immediately before and after the CE phase are ai and at 
respectively. The constant acE is the fraction of the orbital 
energy used to unbind the envelope from the core of the pri- 
mary and A is the ratio between the approximate form of the 
binding energy GMc(Me + Mc)/Ri,l and the true binding 
energy. Throughout our investigation we use A = 0.5, which 
is consistent with the value used by Willems fc Kolb (2004) . 
We take acE ~ 1-0 as our reference model 'A'. 

Once the CE has been ejected from the system, then 
if the binary's stellar components have avoided a collision, 
what remains is a PCEB system. The subsequent evolution 
of the system is driven by orbital AM losses; a combina- 
tion of magnetic braking and gravitational radiation if the 
secondary mass is larger than the mass of a fully convec- 
tive, isolated main sequence star, AfMS,conv, and purely by 
gravitational radiation otherwise. 

The system will become semi-detached once the orbital 
separation has decreased sufficiently so that the secondary 
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Table 1. Model assumptions for acE a^nd magnetic braking. 2 5 



Model 


"CE 






Magnetic braking law 


hA 


1.0 






eqn. ^ 


hCEOl 


0.1 






eqn. ^ 


hCE06 


0.6 






eqn. ^ 


hCE3 


3.0 






eqn. ^ 


hCE5 


5.0 






eqn. ^ 


hPWR05 


eqn. ^ , 


P = 


0.5 


eqn. ^ 


hPWRl 


eqn.(|6ll, 


P = 


1.0 


eqn. ^ 


hPWR2 


eqn.©, 


P = 


2.0 


eqn. (HJ 


rvj2A 
rvj4A 


1.0 
1.0 






eqn. lO, 7 = 2 
eqn. JTl, 7 = 4 



star is brought into contact with its Roche lobe. If the mass 
ratio g<1.3 then mass transfer driven by AM losses will 
ensue. If, however, g>1.3 then the secondary will undergo 
thermal timescale mass transfer (TTMT). Depending on the 
initial mass ratio and evolutionary state, some of these sys- 
tems may reappear as AM-driven CVs. It is thought that as 
many as 40 per cent (|Kolb fc Willemj|2005^ of the present 
population of CVs may have formed in this way. The de- 
tection of CNO abundance anomalies supports this hypoth- 
esis, placing a lower limit of 10 per cent on the fraction 
of post-TTMT systems among the observed population of 
CVs (Gansicke et al. 2003). Determination of the parame- 
ters at which the TTMT phase reverts to the standard AM- 
driven phase would require a comprehensive treatment of the 
TTMT phase, and would depend on the fate of the trans- 
ferred material (Schenker 2001). We therefore only consider 
the contribution of purely AM driven CVs in the present 
investigation. 



2.2.1 The Ejection Efficiency Parameter, acE 

In our calculations, we consider the impact of changing the 
common envelope ejection efficiency acE- In addition to con- 
stant, global values of qce we follow Politano & Weiler 
(2007) and consider the possibility that acE could be a func- 
tion of the secondary mass. We adopt 

= (S) ' 

where p is a free parameter. Figure [2] shows how qce varies 
with M2 for p = 0.5, 1.0 and 2.0. The dark-grey shading is 
the range of secondary masses in gPCEBs considered here, 
while the light-grey shading is the range of secondary masses 
found for dCV progenitors. Table [1] summarises our popu- 
lation synthesis models for varying values of acE- In our 
model names, the first portion in lower case denotes the 
form of magnetic braking used. In Table [T]'h' denotes that 
the Hurley et al. (2002) form of magnetic braking is used ac- 
cording to equation ((8]). The label 'CEa;' denotes the value, 
X, of the constant, global value of acE- For example 'CEOl' 
is acE = 0.1. Finally, the acronym 'PWRj/' denotes the use 
of equation (6) with the value, y, of the parameter p, for 
example 'PWR05' is p = 0.5. 




0.5 - 
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Figure 2. Variation of cice with M2 according to eqn. Solid: 
p = 0.5; dashed: p = 1.0; dot-dashed: p = 2.0. Dark-grey shading 
indicates range of secondary masses in gPCEB systems; light-grey 
range of secondary masses for dCV progenitors. 

2.2.2 Magnetic Braking 

We also investigat e how various forms of m agnetic braking 
affect our results. iRappaport et all l|l983l ) give a parame- 
terised form of magnetic braking as 

j„j = -r^,,^ 3.8 X 10-=^° R% s cm' M2 n^ (7) 

where Q, is the spin frequency of the secondary star and 7 
is a free p arameter. If 7 = 4 then ((7| reduces to the form 
derived bv lVerbunt fc Zwaa The dimensionless con- 

stant 77rvj is determined so as to give the correct AM loss rate 
at the upper edge of the period gap, which we determine in 
Section 2.2.3. The prescription given by Hurley et al. (2002) 
is given by 

Jh = -7?h 5.83 X 10-1^ Mq RqI yr (8) 

where Mcnv is the mass of the secondary star's convective 
envelope and R2 is the donor's radius. Thus, magnetic brak- 
ing vanishes for a star with no convective envelope. As with 
equation ((Tjl the value of the constant T^h is determined to 
give the appropriate AM loss rate at the upper edge of the 
period gap. We use equation ([8]) for our reference model. Ta- 
ble 1 summarises our models where we use different forms 
of magnetic braking. The acronym 'rvj2' and 'rvj4' denotes 
that we use the Rappaport et al. (1983) form of magnetic 
braking with 7 = 2 and 4 respectively. 

2.2.3 The Period Gap 

In this section we outline how we determine the position, and 
therefore the width, of the period gap. We then ascertain the 
values of i^rvj and ?7h , which will give the AM loss rate at the 
upper edge of the period gap appropriate for its determined 
width. Mass loss from the donor will cause it to exceed its 
equilibrium radius (see Section 2.2.5) by a factor /conv at 
the upper edge of the period gap. We calculate this value at 
the end of this section. 

We obtain the position of the gap edges by calculat- 
ing a moving average of the cumulative distribution curve 
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shown in Figure [T] Linear fits were performed on tliree seg- 
ments of tfie curve; witliin and on eitfier side of the pe- 
riod gap. Specifically, we consider the curve in the range 
0.22 < logio(Porb/ h) < 0.64. The edges of the period gap 
are where these fits intersect. We find that the upper and 
lower edges are located at 3.0 and 2.0 hours respectively. We 
adopt these values throughout the rest of our investigation. 

The magnetic braking laws described in Section 2.2.2 
are now calibrated so as to give the correct AM loss rate, and 
therefore the correct mass transfer rate, at the upper edge 
of the period gap. Kolb (1996) determined that the mass 
transfer rate at the upper edge of the period gap is ~ 10~^ 
Mq yr~^. For the Rappaport et al. (1983) form of magnetic 
braking we find that rj^vj = 0.28 for 7 = 2. For 7 = 4 we 
find 77rvj ~ 1.5. For the Hurley et al. (2002) prescription, 
J7h = 0.19. 

We now determine the value of /conv. The Roche lobe 
radius of the donor star, R2,l, from Roche geometry 

, 1/3 



R2 



0.44 



/ Ma y 
\M2 + mJ 



M2 + M1. 

(Paczyiiski 1971) and Kepler's law gives 
^1/2 

Porb OC 



M2" 



(9) 



(10) 



Hence the orbital periods of the upper (Pu) and lower {Pi) 
edges of the period gap relate to the corresponding donor's 
radii at these locations, _Rl u and -Rl i as 



Pi 



f Rl.u 
\Rl.i 



3/2 



(11) 



There is no term in M2 as the mass of the secondary star 
remains constant across the period gap. Given that Pu/Pi = 
3/2, then RL,u/Rh,i = ./conv ~ 1-3. 

2.2.4 Mass Transfer m CVs 

The instantaneous mass transfer rate, M2, from the CV 
don or with radi us R2 and Roche lobe radius -Ra.L is given 
by (|Ritteill 19881 ) 



M2 



Mo exp 



^ R2,L-'R2 ^ 



(12) 



where H ~ 10~*-R2 is the photospheric scale height of the 
donor star and A/q « 10~* Mq yr~^. We assume steady 
state mass transfer, i.e. M2 — 0. From equation (|12|l . this 
assumption implies that the donor radius and Roche lobe 
radius move in step 

^(-R2,L -i?2) =0, (13) 

where we have also assumed that R2 — i?2,L <C R2. 

We calculate the mass transfer rate for a given binary 
configuration time step by considering the relative difference 
between the donor radius and Roche lobe radius 5 = {R2 — 
R2,h)/R2 from BiSEPS as a function of trial values of the 
mass transfer rate. We locate the root of the function S{M2) 
by employing a bisection algorithm, and therefore obtaining 
the value of M2 ■ In practice we iterate until we find a value 
of M2 where |5|<10~^. This procedure is similar to that 
used by Kalogera et al. (2003), although our method has 
been developed independently. 



The bisection algorithm encounters numerical problems 
for those CVs that are close to thermal or dynamical in- 
stability, or where the secondary star has evolved signif- 
icantly from the zero-age main sequence (typically where 
M2/M0>1.O). Specifically, the algorithm obtains unphys- 
ically high mass transfer rates at the onset of mass trans- 
fer. Consequently, the evolutionary tracks of these systems 
across the period gap are not calculated. 

We estimate the present day number of dCVs that these 
systems would produce for each model by multiplying the 
mean dCV lifetime with the formation rate of these numer- 
ically unstable CVs (we apply equation (|25p described in 
Section 3.1). The mean dCV lifetime, (rdcv), was calculated 
by evolving a typical CV with a 0.4 Mq donor and a 0.6 Mq 
white dwarf. These typical values were determined by cal- 
culating the probability distribution of all zero-age CVs. We 
find that (rdcv) ~ 1053 Myr. 

The estimated present day number of these dCVs are 
added onto the present day numbers calculated from the 
population synthesis. The fraction of the total dCV model 
population that these numerically unstable CVs contribute 
are as follows. For the Hurley et al. (2002) prescription of 
magnetic braking, the contribution is ^8.0 per cent. For 
the Rappaport et al. (1983) prescription, the contribution 
is < 20.0 per cent when 7 = 2 and < 30.0 per cent when 
7 = 4. 



2.2.5 The Reaction of the Donor to Mass Loss 

We apply a simple analytical scheme that mimics the re- 
action of the donor star to mass loss. The expansion of a 
low-mass donor as a result of mass loss is a well established 
result from full stellar evolution calculations and analyti- 
cal considerations (Stehle et al. 1996). We define the radius 
excess, /, as the factor by which the star has exceeded its 
equilibrium radius, R2,eq, due to mass loss. We parameterise 
the evolution of / according to a power law expression 



f 



R2: 



R2., 



AM^, 



(14) 



where A and P are constants. 

At the upper edge of the period gap, the donor star 
exceeds its equilibrium radius by a factor / = /conv = 1-3. 
At the point where mass transfer just commences at AI2 — 
M2fi, we require that the donor star is in thermal equilib- 
rium i.e. / = /o = 1. Thus we have 



/o = AM^o, 
and 

/conv ~ j4 Af|^ j,j,jj.j, , 

Combining (fT4)) . (fT5)l and (fT6)) yields 



/ = /c. 



M2 



M2. 



where 



ln(/o//c, 



ln(M2,o/M2 



(15) 



(16) 



(17) 



(18) 



The value of the donor's radius calculated by BiSEPS 
is then multiplied by /. We assume furthermore that the 



6 P. J. Davis et al. 



-0.2 - 



Solid: Main sequence 
DoLLcd: 
Short — daslied 
Long-dashed: 
Dot-dashed: 



-0.4 



-0.6 




-0.8 



-0.6 -0.4 
logj„(M,/MQ) 



Figure 3. Mass-radius relations for CVs with different initial 
donor masses. 



that the donors do exceed their equihbrium radius by the 
same amount by the time they become fully convective. As 
such, we ensure that we reproduce the observed width of 
the period gap with a well defined location. As we are in- 
terested in the evolution of systems across the period gap, 
rather than the detailed CV evolution above the gap, we 
believe the pr escription we adopt i s satisfactory. 

Following 'Kin g fc Kolbl (|l995l ) once the donor star de- 
taches from its Roche lobe it re-attains its equilibrium radius 
on the timescale 



1 



/c. 



D 



Ccff — Cad 



(21) 



where C,eS = dlni?2 /dlnM2 is the effective mass-radius index, 
Cad is the adiabatic mass-radius exponent, rj is the angular 
momentum loss rate timescale and the stability denominato r 
D « 1 is given by equation (16) in iKing fc Kolbl l|l995h . 
These constants are evaluated just before the CV becomes 
detached at the upper edge of the period gap. Typically, 
T ~ 10^ yr. In the detached phase the radius excess follows 
the relation 

f(t) = 1 - (1 - /co, 



-t/T 



(22) 



where t is the time since the donor detached from its Roche 
lobe. 



donor's effective temperature, Tofj, is the same as its iso- 
lated, main sequence equivalent (Stehle et al. 1996). The lu- 
minosity, L2 , of the donor star is calculated from the Stefan- 
Boltzmann Law 

La = A-wRi^oT^s- (19) 

where a is Stefan-Boltzmann constant. Ifence L2 will be 
larger by a factor of compared to models without mass 
loss. The mass-radius relation of several CV evolutionary se- 
quences calculated with this prescription is shown in Figure 

El 

In order for the lower edge of the period gap to be in the 
correct location of Pi — 2.0 h, we disrupt magnetic braking 
when M2 = Af2,conv = 0.17 Mq. 

Based on homology relations, Stehle et al. (1996) de- 
rived a first order differential equation that describes how 
the donor star's radius reacts to mass loss. This scheme 
shows that CVs with different initial system parameters at 
the onset of mass transfer quickly converge to a single, uni- 
form evolutionary track, consistent with the well defined po- 
sition of the period gap. 

Our initial scheme to describe the reaction of the donor 
to mass loss in terms of the mass loss timescale Tj^j^ and its 
Kelvin-Helmholtz timescale tkh, was closely based on that 
of Stehle et al. (1996), according to 

/ = + 1, (20) 

where the value of the constant k can be determined by stip- 
ulating that / = /conv = 1-3 at the upper edge of the period 
gap. However, we encountered numerical instabilities due to 
feedback between the donor radius and the mass transfer 
rate. The scheme we use instead does not reproduce conver- 
gence of CV evolution onto a single track, but does ensure 



3 RESULTS AND ANALYSIS 

We now present the results of our population synthesis cal- 
culations. These are summarised in Table [2] which lists the 
present day population of dCVs, gPCEBs and the ratio 
dCV:gPCEB in the period gap. For completeness we also list 
the total present day formation rates of all PCEBs, and CVs 
above and below the period gap. In all our models, we also 
consider different forms of the IMRD, with — —0.99, 0.0 
or 1.0. We also determine the secondary mass from the same 
IMF as the primary star. Our reference IMRD is n{q) = 1. 

The excess of dCVs over gPCEBs is illustrated in the 
top panel of Figure a prominent peak within the period 
gap in the orbital period distributions of the combined dCV 
and gPCEB populations. These have been calculated from 
our reference model hA, for different assumptions on the ini- 
tial distribution of the secondary mass. The cut-off in the 
number of systems at 2.0 hours is because the smallest sec- 
ondary mass considered in both the dCV and gPCEB sys- 
tems is 0.17 Mq, as such systems become semi-detached at 
2.0 hours. The bottom panel of Figure [4] shows the sepa- 
rate orbital period distributions of dCVs (red) and gPCEBs 
(blue) for our reference model hA, and for n(g) = 1. The 
ratio dCV:gPCEB within the period gap gives an indication 
of the height of the peak. 



3.1 Varying qce: Models hA to hCE5 

We begin by examining how the present day populations of 
dCVs and gPCEBs within the period gap are affected by 
varying the value of the global parameter acE for n{q) = 1 
(Figure [5]). An understanding of these trends will yield an 
understanding of the overall trend of dCV:gPCEB with acE- 
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Figure 4. Top panel: orbital period distribution of dCVs and 
gPCEBs combined for our reference model hA, and for dilferent 
initial mass ratio distributions. Long-dashed: n(q) oc q; short- 
dashed: n{q) = 1 (our reference IMRD); dot-dashed: where M2 
is picked independently from the same IMF as Mi according to 
equation |[TJ; dotted: n{q) oc q—0 S9 (inset). The peak apparent 
in all the distributions is due to the population of dCVs; the 
gPCEB orbital period distributions are increasing monotonically 
with Porb- Bottom panel: orbital period distribution of dCVs 
(red) and gPCEBs (blue) for our reference model hA, and for 
n(g) = 1. 



3.1.1 The gPCEB Population 

The left panel of Figure [5] shows that increasing the value 
of QcE from 0.1 to 1.0 increases the present day number of 
gPCEBs within the period gap (from 1.7 x lO'^ to 1.0 x 10*^ 
for our reference IMRD). This levels off (here to 1.1 x lO'') 
for acE > 1.0. 



To understand this consider Figure [G] where the dashed 
box CDEF defined by 0.17 < A'h/MQ < 0.36 and 
logiQ (Pf ) < Porb/h < logj^Q (Pu) shows the location of 
gPCEBs in the period gap. We designate this region as K. 
Not only will 5R consist of gPCEBs that formed there, but 
also of gPCEBs that evolve into K from longer orbital pe- 
riods. If tu ~ 10 Gyr is the maximum time allowed for the 
gPCEB to reach the upper edge of the period gap and tcR is 
the AM loss timescale, J/J, due to gravitational radiation, 
then the largest orbital period from which a gPCEB can 
evolve into the period gap within the lifetime of the Galaxy, 
PjJ", can be solved from 



tGR(PsJ) - tGR(Pu) 



(23) 



(e.g. Kolb & Stehle 1996). The location of is also shown 
in Figure [S] Therefore, the total area from which gPCEBs 
can evolve into K within the Galactic lifetime is given by 
the region ABCD. 

The whole population of gPCEBs is also bound within 
an upper limit, Pj^ceb (shown in Figure [T]), and the Roche 
lobe-filling limit (RLFL) which is the orbital period where 
systems become semi-detached. As the value of acE is de- 
creased the boundary Pj^ceb ^^'^ population bound 
within will shift to shorter orbital periods. This is because, 
for low ejection efficiencies, more orbital energy is required 
to eject the envelope from the system. The resulting PCEB 
(and therefore gPCEB) population will lie at shorter orbital 
periods. 

We find that the fiux of gPCEBs from ABCD into K 
is larger than the birthrate of gPCEBs within K by a fac- 
tor of ~ 7. If the typical lifetime of a gPCEB within 5R is 
(-TgPCEB) ~ 380 Myr and the fiux of gPCEBs from ABCD 
into K is Psr then the present day number of gPCEBs within 
the period gap, A^gPCEB, can be approximated as 



Fm X (Tg: 



(24) 



Clearly Pj} will depend on the formation rate of gPCEBs in 
ABCD. Figure [7] shows the formation rate of gPCEBs as a 
function of logjQ (Porb/d) for n{q) — 1, and for qce = 0.1 
(red), 1.0 (black) and 5.0 (green). The distribution function 
is fiat-topped and tails off towards longer orbital periods. As 
Qce increases and the whole population of gPCEBs shifts 
to longer orbital periods, the formation rate of the gPCEBs 
between logio(Pu) and Pj^ (inset)- and therefore in ABCD- 
changes as a consequence; at its smallest for qce =0.1 but 
with little change for qce ^1-0. 

From equation (I24p the trend in A^'gpcEB with qce will 
therefore mirror that of the birthrate of gPCEBs in ABCD. 
This is indeed the case as shown in the left panel of Figure 

El 



3.1.2 The dCV Population 

As with the present day population of gPCEBs, there is 
also an increase in the present day number of dCVs associ- 
ated with an increase in qce from 0.1 to 1.0, as shown in 
the middle panel of Figure [5] For our reference IMRD, the 
number of dCVs increases from 2.2 x 10*^ to 3.5 x 10^ In 
contrast to the gPCEB population however, the number of 
dCVs decreases with a further increase in the value of qce- 
If the formation rate of dCVs is Pdcv and their average 
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Table 2. Present day population of dCVs, gPCEBs and the ratio dCVigPCEB in the period gap for various values of «CE a-nd various 
forms of magnetic braking. Also quoted is the present day formation rate of all PCEBs, and of CVs above and below the upper edge 
of the period gap. 



Model 


Number of 


Number of 


Formation rate of 


Formation rate 


Formation rate 


dCV:gPCEB 




dCVs 


gPCEBs 


PCEBs (yr-i) 


of CVs > 3.0 h 


of CVs < 3.0 h 












(yr-i) 


(yr-i) 




n{q) oc < g < 1 


hCEOl 


5.7 X 10'* 


8.7 X 10^ 


2.1 X 10"* 


6.7 X 10-5 


3.7 X 10-5 


6.6 


hCE06 


1.5 X 10^ 


6.1 X 10* 


1.5 X 10-3 


1.8 X 10-* 


2.3 X 10-* 


2.5 


hA 


1.5 X 10^ 


7.4 X 10* 


2.0 X 10-3 


1.8 X 10"* 


2.7 X 10-* 


2.0 


hCE3 


9.5 X 10* 


8.2 X 10* 


3.1 X 10-3 


1.1 X 10-* 


2.4 X 10-* 


1.2 


hCE5 


5.9 X 10" 


7.9 X 10* 


3.5 X 10-3 


7.7 X 10-5 


2.0 X 10-* 


0.7 


hPWR05 


1.6 X 10^ 


5.2 X 10* 


1.6 X 10^3 


1.8 X 10-* 


1.9 X 10-* 


3.1 


hPWRl 


1.4 X 10^ 


2.5 X 10* 


1.3 X 10-3 


1.7 X 10-* 


1.0 X 10-* 


5.6 


hPWR2 


1.0 X 10^ 


3.3 X 103 


1.1 X 10-3 


1.4 X 10-* 


2.2 X 10-5 


30.3 


rvj2A 


2.2 X 10^ 


7.4 X 10* 


2.0 X 10-3 


2.5 X 10-* 


2.8 X 10-* 


3.0 


rvj4A 


2.4 X 10^ 


7.4 X 10* 


2.0 X 10-3 


3.2 X 10-* 


2.8 X 10-* 


3.2 


n(q) = 1.0, < g < 1 (our reference IMRD) 


hCEOl 


2.2 X lO'' 


1.7 X 10^ 


9.8 X 10-3 


3.1 X 10-3 


7.0 X 10-* 


12.9 


hCE06 


3.8 X lO*' 


8.6 X 10^ 


5.5 X 10-2 


5.6 X 10-3 


3.3 X 10-3 


4.4 


hA 


3.5 X lO*' 


1.0 X 10® 


7.2 X 10-2 


5.0 X 10-3 


3.3 X 10-3 


3.5 


hCE3 


2.0 X lO'' 


1.1 X 10® 


1.0 X 10-1 


2.7 X 10-3 


2.5 X 10-3 


1.8 


hCE5 


1.0 X 10'' 


1.1 X 10® 


1.2 X 10-1 


1.5 X 10-3 


1.9 X 10-3 


0.9 


hPWR05 


4.1 X lO*' 


7.7 X 10^ 


6.8 X 10-2 


5.2 X 10-3 


2.8 X 10-3 


5.3 


hPWRl 


3.9 X lO*' 


4.4 X 10^ 


6.5 X 10-2 


5.2 X 10-3 


1.9 X 10-3 


8.9 


hPWR2 


3.1 X 10** 


7.2 X 10* 


6.1 X 10-2 


4.8 X 10-3 


6.8 X 10-* 


43.1 


rvj2A 


5.5 X 10" 


1.0 X 10® 


7.2 X 10^2 


7.6 X 10-3 


3.6 X 10-3 


5.5 


rvj4A 


6.0 X 10" 


1.0 X 10® 


7.2 X 10^2 


1.0 X 10-2 


3.9 X 10-3 


6.0 


n{q) oc q, < q < 1 


hCEOl 


1.8 X lO'' 


6.4 X 10* 


1.1 X 10-2 


3.2 X 10-3 


3.1 X 10-* 


28.1 


hCE06 


2.3 X lO*' 


2.7 X 10^ 


5.8 X 10-2 


4.2 X 10-3 


1.4 X 10-3 


8.5 


hA 


2.0 X 10*^ 


3.1 X 10^ 


7.5 X 10-2 


3.6 X 10-3 


1.1 X 10-3 


6.5 


hCE3 


9.9 X 10^ 


3.6 X 10^ 


1.1 X 10-1 


1.7 X 10-3 


6.1 X 10-* 


2.8 


hCE5 


4.1 X 10^ 


3.7 X 10^ 


1.2 X 10-1 


7.0 X 10-* 


4.5 X 10-* 


1.1 


hPWR05 


2.3 X 10® 


2.5 X 10^ 


7.7 X 10-2 


3.7 X 10-3 


1.1 X 10-3 


9.2 


hPWRl 


2.3 X 10® 


1.6 X 10^ 


7.8 X 10-2 


3.7 X 10-3 


9.1 X 10-* 


14.4 


hPWR2 


2.0 X 10® 


3.2 X 10* 


7.8 X 10-2 


3.5 X 10-3 


5.3 X 10-* 


62.5 


rvj2A 


3.2 X 10® 


3.1 X 10^ 


7.5 X 10-2 


5.7 X 10-3 


1.5 X 10-3 


10.3 


rvj4A 


3.6 X 10® 


3.1 X 10^ 


7.5 X 10-2 


8.9 X 10-3 


1.9 X 10-3 


11.6 


M2 from IMF according to cqn. JTJ 


hCEOl 


1.8 X 10® 


4.0 X 10^ 


6.5 X 10-3 


1.9 X 10-3 


1.7 X 10-3 


4.5 


hCE06 


4.7 X 10® 


2.9 X 10® 


4.7 X 10-2 


4.8 X 10-3 


1.1 X 10-2 


1.6 


hA 


5.0 X 10® 


3.5 X 10® 


6.5 X 10-2 


4.8 X 10-3 


1.3 X 10-2 


1.4 


hCES 


3.1 X 10® 


3.9 X 10® 


1.0 X 10-1 


3.1 X 10-3 


1.2 X 10-2 


0.8 


hCE5 


1.8 X 10® 


3.7 X 10® 


1.2 X 10-1 


2.1 X 10-3 


1.0 X 10-2 


0.5 


hPWR05 


5.2 X 10® 


2.4 X 10® 


4.4 X 10-2 


4.9 X 10-3 


9.0 X 10-3 


2.2 


hPWRl 


4.7 X 10® 


1.2 X 10® 


2.9 X 10-2 


4.6 X 10-3 


4.3 X 10-3 


3.9 


hPWR2 


3.1 X 10® 


1.4 X 10^ 


1.7 X 10-2 


3.6 X 10-3 


8.2 X 10-* 


22.1 


rvj2A 


6.2 X 10® 


3.5 X 10® 


6.5 X 10-2 


6.3 X 10-3 


1.4 X 10-2 


1.8 


rvj4A 


6.8 X 10® 


3.5 X 10® 


6.5 X 10^2 


7.6 X 10-3 


1.4 X 10-2 


1.9 
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Figure 5. From left to right: the present day population of gPCEBs in the period gap (blue) and dCVs (red) as a function of ocEi 
and the ratio dCVigPCEB as a function of acE- Long-dashed: n{q) oc q; short-dashed: n{q) = 1 (our reference IMRD); dot-dashed: 
where M2 is picked independently from the same IMF as A/i according to equation JTJ; dotted: n{q) oc q—0 9S (insets). 



lifetime is (rdcv), then the present-day population of dCVs 
can be approximated by 

A^dcv ~ Bdcv X (Tdcv)- (25) 

If we now assume that the populations above the period 
gap of pre-CVs (with M2 > MMS,conv) and CVs are in a 
steady state, we then have for their respective formation 
rates Bprocv and Bcv 

Bprccv ~ Bcv ~ Bdcv- (26) 

The formation rate of dCVs is therefore linked to that of 
these pre-CVs above the period gap. 

Figure |8] shows the present-day formation rate of pre- 
CVs for n{q) = 1 on the M2 — logio(Porb) plane, for 
acE = 0.1 (top panel), 1.0 (middle panel) and 5.0 (bottom 
panel), where the colour bar at the top indicates the forma- 
tion rate in Myr~^ per bin area. The lower boundary in each 
of the populations is the RLFL, while the upper boundary is 
given by those PCEBs where the time taken for the primary 
progenitor to fill its Roche lobe plus the time taken for the 
secondary to subsequently come into contact with its Roche 
lobe is 10 Gyr (i.e. the Galactic lifetime). The shape of the 
upper boundary is a consequence of the AM loss mechanism; 
gravitational radiation for systems with AI2 < Mms.cohv or 
a combination of magnetic braking and gravitational radia- 
tion for systems with M2 > Mms.cohv 

Inspecting the upper and middle panels of Figure [8] 
shows that there is an increase in the formation rate of pre- 



CVs from QfCE =0.1 to QcE = 1-0. This is a consequence of 
the whole population of PCEBs being shifted to longer or- 
bital periods for larger ejection efficiencies, and hence more 
systems surviving the CE phase. For qce = 5.0 in the lower 
panel, however, there is an overall decrease in the formation 
rate of pre-CVs. Progenitor systems with even shorter ini- 
tial orbital periods than those for models hCEOl or hA can 
now survive the CE phase. Such progenitors result in pre- 
CVs that formed through a case B CE phase with a naked 
helium star remnant. These systems lie below the solid line 
in the bottom panel of Figure [S] while those above the line 
are those pre-CVs that formed through a case C CE phase. 
The latter systems have ended up at much longer orbital 
periods than for models hCEOl and hA. The solid line rep- 
resents those pre-CVs whose progenitor primaries just filled 
their Roche lobes on the base of the AGB. It is the disap- 
pearance of case C remnants and the appearance of the less 
abundant case B remnants that causes the decrease in the 
total pre-CV formation rate. 

From equations (|25p to (|26|) this trend in Bprocv should 
be mirrored in the present day number of dCVs, which is the 
case as shown in the middle panel of Figure [5] 

3.1.3 Overall Trends 

These trends in the populations of gPCEBs and dCVs are 
the same for all forms of the IMRD function, and where the 



10 p. J. Davis et al. 



-0.6 




Figure 6. The dashed box, CDEF, shows the region 5R defined 
by 0.17 < Afa/Mo < 0.36 and logio(P.c) < logio(Porb/d) < 
logio(Pu). The gPCEBs within 3? will fill their Roche lobes within 
the period gap. The largest orbital period a gPCEB can have that 
will evolve into the period gap within the lifetime of the Galaxy 
is given by the boundary . Thus the total area from which 
gPCEBs can flow into the period gap within the Galactic lifetime 
is the region ABCD. 
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Figure 7. The formation rate of gPCEBs as a function of 
l°giO (^orb/d) for n{q) = 1 and for acE = 0.1 (red), 1.0 (black) 
and 5.0 (green). Shown in the inset are the period gap boundaries 
\og\oPi (which is also the RLFL for 0.17 secondaries) and 
logioPu, and PjJ". The upper boundary P^q^s '^CE = 5.0 is 
shown in the main plot as the dashed line. Note that P^ipppj, and 



in fact the whole distribution function, 
as acE decreases. 



shifts to shorter periods 
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Figure 8. The formation rate of pre-CVs for n{q) = 1 on the 
A/2 — logioC^orb) plane for acE = 0.1 (top panel), 1.0 (middle 
panel) and 5.0 (bottom panel). Pre-CVs that lie above the solid 
line in the bottom panel formed through a case C CE phase, while 
those that lie below it formed through a case B CE phase with a 
naked helium star remnant. The colour bar at the top indicates 
the formation rate (Myr^^) per bin area dAfa dlog]^Q(Pojb)- 



secondary mass is determined independently from the IMF 
according to IT}. The largest population of gPCEBs and 
dCVs is obtained when the secondary mass is determined 
for this latter case. The smallest population numbers are 
found for n{q) oc q~'^^^^ . 

Combining the trends of dCVs and gPCEBs gives a de- 
crease in the ratio of dCVigPCEB with increasing acE, seen 
in the right-most panel of Figure [S] For our reference IMRD 
function, dCV:gPCEB decreases from 12.9 for qce = 0.1 to 
0.9 for Qce ~ 5.0. Again this trend is the same for all forms 
of the IMRD function. We note that the ratio dCVigPCEB 
decreases with decreasing value of u. To understand why, 
note that the dCVs will have formed from CVs filling their 
Roche lobes at long orbital periods above the period gap. 
This is in contrast to gPCEBs in the period gap which 
have secondary masses in the range 0.17 < M2/M0 < 0.36. 
Thus, for a given primary mass, the dCV progenitor sys- 
tems will have, on average, more massive secondaries than 
found in gPCEBs. The IMRD n{q) oc q will mostly favour 
the formation of systems with more massive secondaries, 
thus mostly favouring the formation of dCV progenitors over 
the gPCEB progenitors. Thus we obtain the largest values 
of dCV:gPCEB. This is in contrast with the IMRD of the 
form n{q) oc g""'^^ which most favour the formation of sys- 
tems with less massive secondaries, thereby mostly favour- 
ing the formation of gPCEB progenitors over dCV progeni- 
tors. We therefore obtain the smallest values of dCV:gPCEB 
out of the three forms of the IMRD. The smallest values of 
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dCV:gPCEB are obtained when the secondary mass is de- 
termined independently from the primary mass. 

3.2 acE as a function of Secondary Mass: models 
hPWR05 to hPWR2 

In models where qce is determined from equation (|6]) both 
the present day populations of gPCEBs and dCVs decrease 
as the power index p is increased. For our reference IMRD 
function, the left panel of Figure [9] shows that the number 
of gPCEBs decreases from 7.7 x lO'^ to 7.2 x 10*. The middle 
panel shows that the number of dCVs decreases slightly from 
4.1 X 10*^ to 3.1 X lO''. As with the population of gPCEBs 
for models hA to hCE5, this trend depends on the flow rate 
of gPCEBs into the region K, and how this changes as the 
whole population of PCEBs is shifted with acE- 

3.2.1 The gPCEB Population 

The dark-grey area in Figure [2] indicates the range of sec- 
ondary masses for the gPCEB systems. As p increases from 
0.5 to 2, the average corresponding value of acE for these 
systems decreases from ~ 0.5 to ~ 0.1. As a consequence 
the population of gPCEBs will shift to shorter orbital peri- 
ods with increasing p. This in turn will change the formation 
rate of gPCEBs in the region ABCD and hence the flux Fu. 

As can be seen in the left panel of Figure O a value of 
acE decreasing from ~ 0.5 to ~ 0.1 leads to the decrease in 
the present-day number of gPCEBs in the period gap. This 
corresponds to the decrease in the present-day number of 
gPCEBs with increasing p seen in the left panel of Figure |9l 

3.2.2 The dCV Population 

As with the dCVs for models hA to hCE5, the population 
of dCVs for hPL05 to hPL2 depends on the formation rate 
of pre-CVs with M2 > AfMS,conv above the period gap. The 
light-grey region in Figure [5] shows the range of secondary 
masses for dCV progenitor systems. For this range in sec- 
ondary mass there is, on average, little change in qce. In- 
deed, the dominant secondary mass group for pre-CV pro- 
genitors is A^2 ^0.6 M0; here an increase in p from 0.5 to 
2 corresponds to a decrease in acE from ~ 0.7 to ~ 0.3. 
Hence from the middle panel of Figure [S] a slight decrease 
in Ndcv can be expected. This is indeed the case as shown 
in the middle panel of Figure (5] 

3.2.3 Overall Trends 

The trend in the present day populations of gPCEBs and 
dCVs with p is the same for all forms of the IMRD function, 
and where the secondary mass is determined from the same 
IMF as the primary mass. As before, the largest population 
of dCVs and gPCEBs occurs when the secondary mass is 
determined from the same IMF as the primary. The smallest 
populations occur when n(g) oc 

As a consequence of the behaviour discussed in Sec- 
tions 3.2.1 and 3.2.2, the ratio dCV:gPCEB increases as p 
increases. For our reference model, dCV:gPCEB increases 
from 5.3 for p = 0.5 to 43.1 for p = 2. The variation of 



dCVigPCEB with p is shown in the most right hand panel 
in Figure [9] 

Note that these ratios are very large for p = 2 as PCEBs 
with low-mass donors will encounter a merger. Hence the 
vast majority of all CVs will form above the period gap, 
with very few gPCEBs within the period gap. Politano & 
Weiler (2007) suggested that such a strong dependence of 
Qce on M2 could explain the lack of CVs with sub-stellar 
secondaries. However, we note that a lack of calculated 
gPCEBs with low-mass secondaries appears to be inconsis- 
tent with observations (Ritter & Kolb 2003, catalogue Edi- 
tion 7.7 (2006)), which show a modest number of PCEBs 
with M2 « 0.17 Mq. We therefore believe that the p = 2 
model is not a realistic one. As with models hA to hCE5, 
dCV:gPCEB decreases with decreasing values of v. The rea- 
sons for this are the same as those described in Section 3.1. 
The smallest value of dCV:gPCEB occurs when the sec- 
ondary mass is determined independently from the primary 
mass. 



3.3 Alternative Magnetic Braking Laws: Models 
rvj2A and rvj4A 

We now consider the impact of different forms of magnetic 
braking on gPCEBs and dCVs. Increasing 7 from 2 to 4 
increases the present day number of dCVs from 5.5 x lO'' 
for 7 = 2 to 6.0 X 10® for 7 = 4, for our reference IMRD. 
The Rappaport et al. (1983) prescription of magnetic brak- 
ing gives a larger present day population of dCVs than the 
Hurley et al. (2002) prescription, between any IMRD func- 
tion; for our reference model hA, with n(g) — 1, the present 
day number of dCVs is 3.5 x 10®. 

The pre-CV progenitor systems of CVs that form above 
the period gap will initially be driven by magnetic brak- 
ing. Increasing the strength of magnetic braking will shift 
the upper limit Pp^^cy to longer orbital periods for systems 
with Af2 > Afivis.conv Thus more systems will become semi- 
detached within the lifetime of the Galaxy. From equations 
(|25|) to (|26|) . this will result in an increase in the formation 
rate of pre-CVs above the period gap, and therefore increase 
the formation rate of dCVs. In turn, this will increase the 
present day population of dCVs. 

Note that the number of present day gPCEBs does not 
change with magnetic braking law. In our model, the mass 
at which isolated stars become fully convective is 0.35 Mq 
(Hurley et al. 2002). Thus, the evolution of gPCEBs with 
Af2 < 0.35 M0 will be driven by gravitational radiation only. 
As we are only considering PCEBs with secondaries in the 
range 0.17 < M2/MQ < 0.36, only very few systems will 
evolve via magnetic braking while the larger majority will 
be driven via gravitational radiation. 

These trends in the population of dCVs and gPCEBs 
result in an increase in dCV:gPCEB with increasing mag- 
netic braking strength. For our reference IMRD function, 
and where qce ~ 1-0, the ratio increases from 3.5 for the 
Hurley et al. (2002) form of magnetic braking, to 6.0 for 
the Rappaport et al. (1983) prescription, 7 = 4. The largest 
values of dCV:gPCEB occur when the IMRD function is 
71(g) (X q, while the smallest occur when n{q) oc g~°'^^. 
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Figure 9. From left to right: the present day population of gPCEBs in the period gap (blue) and dCVs (red) as a function of the 
power index p, and the ratio dCVigPCEB as a function of p. Long-dashed: n{q) oc q; short-dashed: n{q) = 1; dot-dashed: where M2 is 
picked independently from the same IMF as Mi according to equation JTJ; dotted: n{q) oc q—0-99 (insets). 



4 DISCUSSION 

4.1 Observational Predictions 

In this section we consider the number of present day dCVs 
and gPCEBs, and hence the ratio dCV:gPCEB, given the 
commonly assumed form of the IMRD and values of acE- 
We can then gauge whether such an excess of dCVs over 
gPCEBs can be determined observationally. 

Numerical calculations of the CE phase suggest that 
acE = 0.13 to 0.6 (Iben fc Livid 1 19931 ). while the IMRD 
seems to be flat in q, i.e. n{q) = 1 (Goldberg, Mazeh & 
Latham 2003; Mazeh et al. 1992; Duquennoy & Mayor 1991). 
From Table [5] we can estimate that the number of dCVs 
currently occupying the Galaxy is (2.2 — 3.8) x 10^, while 
the number of gPCEBs is in the range (1.7—8.6) x 10^, giving 
a factor of ~ 4 to ~ 13 more dCVs than gPCEBs. Hence we 
predict a prominent peak in the orbital period distribution 
of short orbital period WDMS systems as shown in Figure 
1101 In determining these predictions, we have assumed that 
magnetic braking takes the form shown in equation (|8]). 

Of course, the AM loss rate associated with magnetic 
braking is not very well constrained, and indeed it has 
been suggested that the magnetic braking strength is ~ 100 
weaker than the forms we have used in this investigation 
(Andronov et al. 2003). However, we have not considered 
the forms of magnetic braking that include a saturation 
point if the secondary star has a spin frequency larger than 



some critical value ficrit (Sills et al. 2000; Ivanova & Taam 
2003). This is because, in the case of the Sills et al. (2000) 
form of magnetic braking, the AM loss rate would be too 
small to sufficiently drive the donor star out of thermal 
equilibrium, and consequently produce a period gap. Fur- 
thermore, the Sills et al. (2000) form of magnetic braking 
predict mass transfer rates at the upper edge of the period 
gap that are too low compared to observations, and as such 
would not explain the presence of nova-like CVs above the 
period gap (e.g. Kolb 2002). The Ivanova & Taam (2003) 
form of magnetic braking on the other hand, which is ~ 2 
to 3 times weaker than the Verbunt & Zwaan (1982) form, 
would drive the donor star out of thermal equilibrium but 
not enough to explain the observed width of the period gap. 

Such forms of magnetic braking are therefore inconsis- 
tent with the disrupted magnetic braking model. Instead, 
alternative mechanisms must exist that combine in strength 
to reach the value that we consider at the upper edge of the 
period gap. 



4.2 The Location of the Period Gap 

We now discuss the impact the location of the period gap 
has on our results. We note that Knigge (2006), using a sim- 
ilar method to ours, determined that the period gap edges 
are 2.15 h and 3.18 h, which agrees with our determined 
period gap width of one hour. We attribute the slight differ- 
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Figure 10. The orbital period distribution of the combined pop- 
ulation of dCVs and gPCEBs for n{q) = 1, acE = 0.6 (dashed 
hne) and «CE = 0.1 (dotted hne). 



ence between the two locations of the period gap to the two 
different catalogues used. Knigge (2006) used the Ritter & 
Kolb (2003) catalogue (Edition 7.6), while we use the latest 
Edition 7.7. 

If the period gap was located according to Knigge 
(2006), this would have the effect of increasing the popu- 
lation of gPCEBs within the period gap, and consequently 
decrease the ratio dCV:gPCEB. The higher values of Pu 
and Pe would shift the mass range of gPCEBs to 0.19 < 
M2/M0 < 0.40, and therefore a larger fraction of gPCEBs 
(those with 0.35 < M2/M0 < 0.40) would be driven by 
magnetic braking. In turn, the flow rate of gPCEBs from 
longer orbital periods into J? would increase. We also point 
out that a similar uncertainty is introduced with the value of 
MMS.conv, as this is not well constrained by stellar evolution 
models. 

In their studies on circumbinary disk driven CVs 
Willems et al. (2005) considered the CV population below 
the gap, while Willems et al. (2007) studied the population 
above the gap. These authors never attempted to systemati- 
cally model, or to achieve a best fit model for, the combined 
population. This is why their adopted period gap boundaries 
of 2.25 and 2.75 hours, respectively, result in an implausibly 
narrow gap when taken together. 

Nonetheless we investigated how a narrower period gap 
would affect our results by calculating a further population 
model with the gap boundaries set to 2.25 and 2.75 hours. 
We considered acE ~ 1.0 with our standard IMRD, and re- 
calibrated equation (|8]) to give the AM loss rate appropriate 
for the narrower gap width. We find that dCV:gPCEB de- 
creases to 2.1 compared to our model hA value of 3.5. Thus 
the prominence of the peak would decrease somewhat. This 
is because the time for a dCV to cross a narrower gap would 
decrease (we find (t^cv) ~ 438 Myr), decreasing the present 
day population from 3.5 x 10^ for model hA to 1.1 x 10®. 
The number of gPCEBs in the period gap also decreases 



from 1.0 X 10" for model hA to 5.2 x 10^ because the sys- 
tems occupy a smaller area in period space. 

We emphasise that a gap width of half an hour is a 
somewhat extreme assumption. Despite this we still obtain 
a clear excess of dCVs over gPCEBs in the period gap, and 
hence a prominent peak would still be detected. Indeed, a 
prominent peak would still be detected for any reasonable 
gap width. 



4.3 TTMT Systems 

Observational investigations have shown that in many CVs 
above the period gap the donor is unlikely to be an un- 
evolved main sequence donor (e.g. Beuermann et al. 1998). 
Instead, the spectral types of the donor are later than those 
of an isolated star of the same mass, while donors at orbital 
periods longer than 5 to 6 hours are nuclear evolved (Patter- 
son et al. 2005; Knigge 2006). Theoretical calculations also 
predict this effect (Baraffe & Kolb 2000). It has been sug- 
gested that such systems had an initial TTMT phase where 
the secondary star was much more massive than the white 
dwarf (Schenker & King 2002). Kolb & Willems (2005) find 
that as much as 40 per cent of zero-age CVs form with a 
donor that has evolved more than half way through its main 
sequence lifetime. Gansicke et al. (2003) have detected CNO 
abundance anomalies consistent with nuclear evolved donors 
among ~ 10 — 15 per cent CVs in a sample observed in the 
ultraviolet, corroborating the hypothesis that a significant 
fraction of the known CV population may be post-TTMT 
systems. 

Such CVs, if they appear after their initial TTMT phase 
above the period gap, may contribute to the present day 
population of dCVs. Thus our values of dCV:gPCEB, which 
are calculated from the evolution of purely AM driven CVs, 
are lower limits. To calculate the contribution of post-TTMT 
CVs to the total dCV population would require a compre- 
hensive treatment of the TTMT phase including the fate 
of the transferred material, which determines the point of 
re-appearance of the system as a 'normal' AM driven CV. 



4.4 Magnetic CVs 

We now consider how the contribution of magnetic CVs 
(mCVs) may affect the ratio dCV:gPCEB within the pe- 
riod gap. Intermediate polars (IPs) are thought to evolve 
as non-magnetic CVs and hence become detached at the 
upper edge of the period gap due to the disruption of mag- 
netic braking. Their contribution to the population of dCVs 
is therefore already taken into account in our calculations, 
as are their progenitors to the gPCEB population. 

Polars, on the other hand, may not become detached 
at the upper edge of the period gap, as it has been sug- 
gested that magnetic braking for these systems is estimated 
to be 1 to 2 orders of magnitude weaker than for non- 
magnetic CVs (Li, Wu & Wickramasinghe 1994; Li, Wu & 
Wickramasinghe 1994b). This is supported observationally 
from their white dwarf temperatures: the accretion-heated 
white dwarfs in polars are consistently colder than those 
in non-magnetic CVs at similar orbital periods, implying 
lower accretion rates, and hence angular momentum loss 
rates (Araujo-Betancor et al. 2005). This difference is es- 
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pecially pronounced above the period gap, where magnetic 
braking is the dominant AM loss agent. 

Furthermore, no progenitor WDMS systems of mCVs 
with a magnetic white dwarf has as yet been detected 
(Liebert et al. 2005). It is possible that such progenitors 
evolve in the same way as progenitors of non-magnetic CVs. 
In such a case, the population of WDMS progenitors of po- 
lars will already be taken into account in our calculations 
and are included in the number of gPCEBs we obtained. 

The overall effect of the above on our results would be to 
reduce the number of dCVs and hence the ratio dCV:gPCEB 
by the fraction of polars among the total CV population. 
From the Ritter & Kolb (2002) catalogue. Edition 7.7 (2006), 
we find that polars contribute ~ 13 per cent to the total CV 
population. 

4.5 PCEB Candidates 

Here we discuss the current situation regarding PCEB 
candidates and the surveys presently searching for them. 
Schreiber & Gansicke (2003) studied 30 well observed 
PCEBs, which were chosen to be representative of PCEBs 
that will begin mass transfer within the Hubble time. These 
systems have orbital periods < 2 days, have mass ratios less 
than 1, and do not include secondaries that are sub-giants. 
They found that the majority of systems contained young 
(<5 X 10* yrs), hot 15 000-22 000 K) white dwarfs. 
This is a consequence of the fact that until recently PCEB 
candidates were almost exclusively identified in blue colour 
surveys such as the Palomar Green survey. Thus, systems 
containing cool white dwarfs and/or early type companion 
stars have been missed in most previous work, with the ex- 
ception of a few systems identified as nearby large proper 
motion systems (e.g. RRCae) or spectroscopic binaries (e.g. 
V471Tau). 

The potential of observational population studies of 
PCEBs has improved over the past few years dramatically 
through the SDSS. Because of the vast ugriz colour space 
probed by SDSS, coupled with the availability of high- 
quality follow-up spectroscopy of a large number of objects 
with non-stellar colours, has already lead to the identifi- 
cation of more than 1000 WDMS binaries (Silvestri et al. 
2007), most of them being by-products of the search for 
quasars. A complementary program targeting specifically 
WDMS binaries containing cool white dwarfs and/or early 
type companions is underway to compensate the bias against 
such systems in the previous surveys (Schreiber et al. 2007). 
For the first phase, the PCEBs among the full sample of 
WDMS binaries have to be identified through radial ve- 
locity studies. A first effort along these lines has been car- 
ried out by Rebassa-Mansergas et al. (2007) who identified 
18 PCEB candidates from multiple SDSS spectroscopy ob- 
tained for 101 WDMS. Rebassa-Mansergas (2008, in prepa- 
ration) confirmed through follow-up spectroscopy so far six 
of those candidates as PCEBs with orbital periods ranging 
from 164 min to 1048 min. Additional identification work of 
the PCEBs among the SDSS WDMS is currently carried out 
at the VLT (Schreiber et al. 2008) and the WHT (Gansicke 
et al. in preparation). 

It is foreseeable that, sufficient observational effort be- 
ing invested, it will be possible to build up an orbital period 
distribution of close WDMS binaries comprising potentially 



a few hundred systems. While the SDSS PCEB sample will 
not be free of selection effects (see e.g. Rebassa-Mansergas 
et al. 2007; Pretorius, Knigge & Kolb 2006), those biases 
can be modelled to a large degree. It appears hence feasi- 
ble to subject the disrupted magnetic braking model to a 
stringent test by comparing such an observed orbital period 
distribution to the predictions of our work here (Figure [TO}. 



5 CONCLUSIONS 

We have performed population synthesis calculations to ob- 
tain the present day population of two types of white dwarf- 
main sequence star systems within the 2 to 3 hour cata- 
clysmic variable period gap. The first are post-CE binaries 
with secondaries that have masses 0.17 < M2/M0 < 0.36, 
and so will commence mass transfer in the period gap 
(gPCEBs). The second type are systems that were CVs in 
the past, but detached at the upper edge of the period gap 
as a consequence of disrupted magnetic braking, and are 
crossing the period gap via gravitational radiation (dCVs). 

Our calculations were repeated to consider constant, 
global values of the CE ejection efficiency, acE, and cases 
where acE is a function of secondary mass according to 
equation (O . We considered various forms of magnetic brak- 
ing according to equation ((7|) with 7 = 2 and 4, and equation 

We find that there is a prominent peak in the orbital pe- 
riod distribution of the combined dCV and gPCEB popula- 
tion, due to the excess of dCVs over gPCEBs within the pe- 
riod gap. We find that the ratio dCV:gPCEB, which gives an 
indication of the peak's height, decreases with an increasing 
global value of qce , while increasing with increasing value of 
p. These trends are the same for all assumptions on the ini- 
tial secondary mass distribution. The value of dCV:gPCEB 
ranges from 0.5 for model hCE5 where the secondary mass 
is determined independently from the primary mass using 
the same IMF, to 62.5 for model hPWR2 where n{q) oc q. 

We find further that dCV:gPCEB increases with in- 
creasing magnetic braking strength, although only slightly. 
For n[q) = 1, dCV:gPCEB increases from 3.5 for model hA 
(the weakest) to 6.0 for rvj4A (the strongest). 

For all our models we find that the values of 
dCV:gPCEB are largest when the IMRD function has the 
form n(q) oc q i.e. on average systems will have more mas- 
sive secondary stars. The smallest values occur when the 
secondary mass is determined from the same IMF as the 
primary. 

The most likely value of dCV:gPCEB is between ~ 4 
to ~ 13, thus we can expect a significant peak as shown in 
Figure 1 101 We suggest that if such a feature is observed in 
the orbital period distribution of short orbital period WDMS 
binaries, this would strongly corroborate the disruption of 
magnetic braking. 
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